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ABSTRACT 


A  nonlinear  Boussinesq  model  for  the  shoaling  of  ocean  surface  gravity  waves  on 
beaches  is  presented  and  compared  to  second-order  finite  depth  theory.  The  spectral 
Boussinesq  model  of  Freilich  and  Guza  (1984)  for  uni-directional  waves  propagating 
perpendicular  to  a  beach  with  straight  and  parallel  depth  contours  is  extended  to  obliquely 
propagating  waves.  Predictions  of  the  shoaling  evolution  of  a  single  resonant  triad  with 
two  primary  incident  wave  components  driving  a  secondary  wave  component  are 
compared  to  finite  depth  theory  predictions  of  forced  secondary  waves.  Results  for  both 
sum-  and  difference-interactions  are  presented  for  a  range  of  beach  slopes,  incident  wave 
amplitudes,  frequencies  and  propagation  directions. 

The  comparisons  show  that  there  is  a  region  (roughly  between  1 0  and  4  m  depth 
for  typical  swell  amplitudes  and  frequencies)  where  both  theories  predict  very  similar 
growth  of  secondary  wave  components.  Whereas  Boussinesq  theory  typically  predicts 
slightly  smaller  secondary  wave  amplitudes  than  finite  depth  theory,  the  dependence  of 
the  secondary  wave  response  on  the  directional  spreading  angle  of  the  primary  waves 
predicted  by  both  theories  are  in  good  agreement.  However,  pronounced  discrepancies 
between  Boussinesq  and  finite  depth  predictions  are  noted  for  very  low  (infragravity) 
fi-equency  secondary  waves  on  relatively  steep  beaches. 
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I.  INTRODUCTION 


Wind-generated  ocean  surface  gravity  waves  are  the  principal  driving  force  for 
fluid  motions  and  sediment  transport  in  the  littoral  zone.  Accurate  models  for  the 
transformation  of  waves  on  beaches  are  needed  to  predict  alongshore  and  rip  currents, 
undertow,  bedforms  (e.g.,  bars  and  cusps)  and  beach  erosion  and  accretion.  In  general,  as 
waves  shoal,  amplitudes  increase,  wavelengths  decrease  and  propagation  directions 
change  towards  normal  incidence  to  the  beach.  In  addition  to  these  linear  shoaling  and 
refraction  effects,  nonlinear  wave-wave  interactions  cause  transfers  of  energy  to  wave 
components  with  both  higher  and  lower  frequencies. 

Theories  for  nonlinear  wave- wave  interactions  in  deep  water  (kh»l,  where  k  is 
the  wavenumber  and  h  is  the  water  depth)  and  intermediate  depths  (kh=0(l))  are  usually 
based  on  Stokes’  (1847)  perturbation  expansion  for  weak  nonlinearity  (ak«l,  where  a  is 
the  wave  amplitude)  (Phillips,  1960;  Hasselmann,  1962;  Longuet-Higgins  and  Stewart, 
1962).  A  nonlinear  transfer  of  energy  to  higher  frequencies  occurs  through  sum- 
interactions:  two  primary  wind- waves  having  frequencies  co,  and  CO2,  and  wavenumbers 
k,  and  kj,  excite  a  secondary  wave,  with  the  sum-frequency  co,  -t-coj  and  (vector)  sum 
wavenumber  k,  -1-  k^ .  Phillips  (1960)  showed  that  at  this  order  there  are  no  resonances; 
secondary  waves  are  forced  waves  that  do  not  satisfy  the  gravity  wave  dispersion  relation. 
Since  the  amplitudes  are  bounded,  these  forced  waves  are  often  called  “bound  waves”. 

The  predicted  forced  secondary  waves  are  phase  locked  to  the  primary  wind 
waves,  and  typically  distort  a  linear  sinusoidal  wave  profile  to  shapes  with  peaky  crests 
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and  broad  troughs.  Observations  of  high  frequency  pressure  and  velocity  fluctuations  in 
intermediate  depths  agree  well  with  second-order  nonlinear  theory  predictions  (Herbers  et 
al.,  1992;  Herbers  and  Guza,  1991, 1992,  1994). 

In  shallow  water,  the  Stokes-type  perturbation  solutions  are  valid  only  if  the 

3. 

Ursell  number,  Ur  =  TyjT’  small  (Stokes,  1847;  Ursell,  1953).  While  this  assumption 

K  11 

holds  well  for  intermediate  continental  shelf  depths,  it  usually  breaks  down  on  beaches, 
where  Ur  is  0(1).  In  this  shallow  region  (typically  h<5  m),  where  dispersion  is  weak, 
nonlinearities  become  quite  significant,  and  asymmetrical  wave  shapes  develop  with 
steep  forward  faces. 

The  nonlinear  shoaling  region  where  Ur  is  typically  0(1)  is  better  represented  by 

the  Boussinesq  equations  (Boussinesq,  1871),  which  assume  that  nonlinearity  -,  and 

h 

dispersion  (kh)^  are  both  weak  and  of  the  same  order  (i.e.,  Ur=0(l)).  In  Boussinesq 
theory,  continuous  nonlinear  energy  transfers  take  place  between  near-resonant  wave 
triads.  Mei  and  LeMehaute  (1966)  and  Peregrine  (1967)  extended  the  Boussinesq 
equations  to  a  bottom  sloping  gently  in  2  dimensions,  and  these  equations  have  been  used 
extensively  to  model  the  nonlinear  evolution  of  waves  shoaling  on  beaches. 

Freilich  and  Guza  (1984)  developed  a  Boussinesq  model  for  unidirectional  waves 
propagating  over  a  bottom  varying  slowly  in  one  dimension.  The  predicted  shoaling 
evolution  of  wave  frequency  spectra  agrees  well  with  field  observations  (Freilich  and 
Guza,  1984;  Elgar  and  Guza,  1985,  1986;  Elgar  et  al.,  1990).  Liu,  Yoon,  and  Kirby 
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(1985)  developed  a  parabolic  approximation  to  Peregrine’s  equations,  allowing  for  small 
directional  spreading  angles  on  a  two  dimensional  beach,  but  applications  to  a  continuous 
frequency-directional  spectrum  of  waves  have  not  been  reported. 

Presently,  no  uniformly  valid  theory  exists  for  shoaling  waves.  Whereas  fully 
dispersive  finite  depth  solutions  are  singular  at  the  shoreline,  weakly  dispersive 
Boussinesq  solutions  break  down  in  intermediate  depths  (kh=0(l)).  Thus,  both  finite 
depth  and  Boussinesq  solutions  are  needed  to  model  the  shoaling  zone.  For  typical  swell 
amplitudes  (0(1  m))  and  periods  (0(10  s)),  finite  depth  theory  can  be  applied  in  the 
“outer  shoaling  region”,  extending  from  roughly  20-6  m  depth  (Berbers  et  al.,  1992),  and 
Boussinesq  theory  in  the  “inner  shoaling  region”,  between  roughly  7  m  depth  and  the 
outer  edge  of  the  breaker  zone  (e.g.,  Freilich  and  Guza,  1984;  Elgar  and  Guza,  1985). 
Although  formally  the  validity  ranges  of  finite  depth  and  Boussinesq  theories  do  not 
overlap,  for  practical  applications  there  is  a  region  (roughly  5-8  m  depth)  where  both 
theories  are  expected  to  be  reasonably  accurate.  The  errors  in  finite  depth  and  Boussinesq 
model  predictions  in  this  transition  region  are  poorly  understood. 

Difference-interactions  between  two  wave  components  (co,,k,)  and  ((02,k2) 
excite  a  third  wave  at  the  difference-frequency  and  wavenumber,  (o),  -©2»k,  -kj)  (for 
CO,  >  ©2 ),  that  typically  causes  a  “set  down”  of  the  mean  ocean  surface  under  groups  of 
high  waves  (Longuet-Higgins  and  Stewart,  1962).  Difference-interactions  are 
particularly  important  because  they  drive  energetic  nearshore  low  frequency  motions 
(nominally  0.005-0.05  Hz),  known  as  infragravity  waves  or  “surf  beat”  (e.g.,  Munk, 
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1949;  Tucker,  1950;  Herbers  et  al.,  1995a).  Infragravity  waves  cause  seiches  in  small 
harbors  (Okihiro,  et  al.,  1993),  and  are  believed  to  be  the  primary  means  by  which 
sediment  is  suspended  in  the  inner  surf  zone  during  storms  (Holman  et  al.,  1978;  Wright 
et  al.,  1982;  Beach  and  Sternberg,  1988). 

Whereas  incident  wind  waves  and  the  associated  higher  frequency  secondary 
waves  are  usually  dissipated  in  the  surf-zone,  the  low-ffequency,  long-wavelength 
secondary  waves  forced  by  difference-interactions  may  reflect  from  shore.  Longuet- 
Higgins  and  Stewart  (1962)  suggested  that  forced  secondary  waves  at  infragravity 
frequencies  are  released  in  the  surf  zone,  reflect  from  the  beach,  and  subsequently  radiate 
seaward  as  free  waves.  Observations  of  the  directional  properties  of  infragravity  waves 
are  consistent  with  this  mechanism  (Herbers  et  al.,  1995b),  but  quantitative  models  for 
the  nonlinear  transfer  of  energy  to  infragravity  frequencies  are  presently  lacking.  Second 
order  finite  depth  theory  (Hasselmann,  1962,  based  on  the  Stokes  expansion)  predicts 
accurately  the  amplitudes  of  forced  infragravity  waves  in  intermediate  depths  (Herbers  et 
al.,  1994),  but  is  singular  near  the  shoreline,  and  thus  cannot  predict  the  amplitudes  of 
free  infragravity  waves  released  in  the  surf  zone  (Herbers  et  al.,  1995b).  Boussinesq 
theory  may  provide  quantitative  estimates  of  free  infragravity  energy  radiated  from  shore. 

This  paper  extends  the  work  of  Freilich  and  Guza  (1984),  by  examining  both  sum- 
and  difference-interactions  between  obliquely  propagating  shoaling  waves.  The  2- 
dimensional  Boussinesq  equations  for  a  beach  with  straight  and  parallel  depth  contours 
are  derived  in  Chapter  II,  and  solved  for  a  spectrum  of  waves  in  Chapter  III.  The 
predicted  evolution  of  amplitudes  and  phases  for  a  single  wave  triad  are  compared  to 
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predictions  based  on  Hasselmann’s  (1962)  second  order  finite  depth  theory  in  Chapter  IV. 
The  results  are  summarized  in  Chapter  V.  All  Figures  are  contained  in  the  Appendix. 
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II.  GOVERNING  EQUATIONS  AND  APPROXIMATIONS 

In  this  chapter,  Boussinesq  equations  are  derived  for  obliquely  propagating  waves 
on  a  gently  sloping,  impermeable  beach  with  straight  and  parallel  depth  contours. 
Following  Peregrine  (1967),  non-dimensional  variables  will  be  used  as  follows: 


h  = 


hi 

ho 


(x,y,z)=  -^(x*,y*,z*)  t  =  ^  t  * 


(u,v,w)=-^J=:-(u*,v*,w*)  p  =  — ^p* 

Wo  pgho 


where  *  indicates  a  dimensional  variable  and  h^  is  a  representative  water  depth.  The  z 
axis  points  vertically  upwards,  with  z=0  corresponding  to  the  mean  surface.  The  x  axis 
points  onshore,  and  y  points  alongshore  (Figure  1).  The  variables  p,  u,  v,  and  w  denote 
pressure  and  the  velocity  components  in  the  x,  y,  and  z  directions,  respectively.  The  free 
surface  is  defined  by  z  =  Ti(x,t)  where  x  is  the  horizontal  space  vector  (x,y).  The  bottom 

is  defined  by  z  =  -h(x). 

The  equations  of  motion  and  the  continuity  equation  in  non-dimensionalized  form 
(neglecting  viscous  effects)  are: 


— -i-(u-V)u  +  w— -i-Vp=  0  (2-1) 

dt  dz 


1 


(2-2) 


^  /  V7\  ^  ^  ^ 

— +  (u-V)w  +  w —  +  — +  1  =  0 
dz  dz 


_  5w  ^ 
V-U  +  —  =  0 
dz 


id  5  1 

where  V  denotes  the  two  dimensional  gradient  operator  — , —  ,  and  u  indicates  the 

v5x  dyj 

horizontal  velocity  vector  (u,v).  The  surface  and  bottom  boundary  conditions  are  given 


w=— +u  Vti  at  z  =  ri(x,t) 


(2-4a) 


at  z  =  Ti(x,t) 


(2-4b) 


w  =  -u—  at  z  =  -h(x) 
dx  ^  ^ 


(2-4c) 


Irrotational  flow  is  assumed,  which  yields  the  relation: 


=  Vw 


The  continuity  Equation  (2-3)  can  be  vertically  integrated  using  the  kinematic 
boundary  conditions  (2-4a,c)  and  Leibniz’s  Theorem: 


V-  [udz  +  —  =  0 

{  at 
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Following  Peregrine  (1967),  the  variables  t]  ,  p,  and  u  are  expanded  in  terms  of 
the  non-linearity  parameter  8 ,  defined  to  be  the  ratio  of  wave  amplitude  to  water  depth 


T1  =  ST1,  +8^112+... 

p  =  Po +Ep, +8^P2+...  (2-7) 

U  =  8U|  +8^U2+... 

where  Po  =  -z,  is  the  hydrostatic  pressure  of  the  undisturbed  fluid.  In  shallow  water,  the 
vertical  velocity  component  w  is  0(ct)  smaller  than  u,  where  the  dispersion  parameter, 
CT  =  kh ,  is  the  product  of  the  wavenumber  and  the  water  depth; 

w  =  a(8w,  +  8^W2+...)  (2-8) 

Even  waves  with  a  large  oblique  incidence  angle  Oq  relative  to  the  shore-normal 
in  deep  water  will  refract  so  that  the  angle,  6(x),  in  shallow  water  is  small  and  thus 
V  «  u .  Assuming  that  the  waves  incident  from  deep  water  propagate  over  a  seabed  with 
approximately  straight  and  parallel  depth  contours,  the  scaling  of  the  v  component 
follows  from  Snell’s  law: 


where  co  and  k(x)  =(k(x),f)  are  the  non-dimensionalized  wave  frequency  and 
wavenumber  vector.  In  shallow  water  this  relation  reduces  to 
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0(x)w  k(x)h(x)sin0o 

Hence,  the  ratio  between  v  and  u  in  shallow  water  is  0(kh),  and  we  can  write; 


u  =  8U,  +  8  U2+... 


V  =  a(8v,  +  8  V2+. 


(2-10) 


The  independent  variables  are  scaled  accordingly: 


dx  ax'  ’  5y  dy'  '  dz  dz'  ’  di~^  di' 


(2-11) 


Since  the  variables  p.,  u.  and  w.  (i  =  0,1,...)  are  0(l),  the  order  of  the  terms  appears 

explicitly  when  the  relations  (2-7),  (2-8),  (2-10)  and  (2-11)  are  substituted  into  the 
governing  equations. 

The  water  depth  h  is  assumed  to  be  a  slow  function  of  x' . 


h  =  h(x) 


(2-12) 


where  x  =  8x'  is  a  slow  space  variable,  so  that 


a 

ax'  ax 


(2-13) 


Note  that  in  dimensional  coordinates,  the  beach  slope  is  0{az).  Thus,  for  typical 

shoaling  wave  values  of  8  =  0.05,  a  =  0.2,  the  beach  slope  is  assumed  to  be  0(l  0“^ )  . 
This  is  a  realistic  slope  for  a  natural  sandy  beach. 
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First  order  relations  for  the  wave-induced  velocity  and  pressure  fields  are  obtained 


by  substituting  Equations  (2-7),  (2-8),  (2-10)  and  (2-11)  into  Equations  (2-2)  and  (2-5) 


5u, 

dz' 


=  0(8) 


(2- 14a) 


(2- 14b) 


^  =  0(E)  (2-14C) 

It  follows  from  Equation  (2-14)  that  Uj,  Vj,  and  pi  are  independent  of  z' . 
Equation  2- 14c  integrates  to  (using  the  dynamic  free  surface  boundary  condition 
(Equation  (2-4b))): 


p,=ri,  (2-15) 

The  first  order  horizontal  momentum  equations  can  be  obtained  by  substituting 
Equations  (2-7),  (2-8),  (2-10),  (2-11),  and  (2-15)  into  Equation  (2-1).  In  component 
form: 


^1  ,  ^1 
di'  dx' 


=  0(8) 


(2- 16a) 


dt'  dy' 


=  0(8) 


(2- 16b) 
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Similarly,  substituting  Equations  (2-7),  (2-10),  (2,11),  and  (2-13)  in  Equation  (2- 
6)  yields  the  first  order  depth-integrated  continuity  equation: 


^1 

at' 


au, 

+  h-  ' 


ax' 


=  0(£) 


(2-17) 


Substituting  Equations  (2-7),  (2-8),  (2-10),  and  (2-11)  in  Equations  (2-2)  and  (2-5) 
we  obtain  the  second  order  relations 


auj 


aw, 

ax' 


+  0(8) 


(2-1 8a) 


av; 

az' 


aw| 


+  0(8) 


(2- 18b) 


^2 

dz' 


aw. 


+  0(8) 


(2-1 8c) 


Vertically  integrating  the  first  order  continuity  equation  (Equations  (2-3),  (2-10), 


(2-11)): 


ax'  az' 


using  Equations  (2-4c),  and  (2-13)  yields 


w,(z')  =  -(z'  +  h)-^  +  0(8) 
ox 


(2-19) 


Substitution  of  Equation  (2-19)  in  Equation  (2-18a-c)  and  integration  with  respect 
to  z'  (using  the  dynamic  free  surface  boundary  condition  (Equation  2-4b))  yields: 
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(2-20a) 


,  2  ”1  p2 

u,  =  U,(x',y',t')-  z'h  +  — —  - — — i-O(s) 

^  2  J  dx'dy' 


f2 

V2=V,(x',yM')-  z'h  +  —  (2-20b) 

Pi  =n2(>‘'.y’.<')+ z'h+^  (2-20c) 

2  J  dK  dy 

where  Uj  and  Wj  arbitrary  velocity  functions  arising  from  the  integration. 

These  functions  can  be  expressed  in  terms  of  the  depth  averaged  second  order  flows  Uj 
and  v,  by  vertically  integrating  Equations  (2-20a,b): 

1  1 

u,=— ru,dz'=U,(x',y',t'>-h^  ^  +0(8)  (2-21a) 

‘  h+TT^-h  3  dx 


V2=-^r  v2dz'=V2(x',y',t')+^h^  ^  ^  +0(8)  (2-21b) 

h+q"'-*’  3  5x'5y' 

Finally,  substituting  the  perturbation  expansions  (2-7),  (2-8),  (2-10),  and  (2-11), 
the  first  order  relation  (2-15)  and  the  second-order  relations  (2-20a,b,c)  and  (2-21a,b) 
back  into  the  horizontal  momentum  Equation  (2-1)  yields,  componentwise: 


^1  I  I  I  I  u 

ax'  at'  lax'  a'  '  ax'  3  ax'^a' 


0(8^) 


(2-22a) 
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av,  1  2  a'u, 


and  the  corresponding  second-order  depth-integrated  continuity  equation  is  obtained  by 
substituting  Equations  (2-7),  (2-10),  (2-1 1),  and  (2-13)  into  Equation  (2-6): 


X.  5h 


=0(s>)  (2-22C, 
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III.  SPECTRAL  SOLUTIONS 


In  Chapter  II,  an  approximate  set  of  governing  equations  was  derived  (2-22a,b,c) 
for  waves  shoaling  on  a  beach  with  straight  and  parallel  depth  contours.  In  this  chapter, 
these  equations  will  be  solved  for  a  frequency-directional  spectrum  of  waves.  Equations 
(2-22a,b,c)  can  be  expressed  in  terms  of  the  first-  and  (depth-averaged)  second-order 

velocity  potential  functions  (j),  and  ^2  • 


u,= 


5x' 


;  v,= 


dy' 


U2=-^+0(£);  V2=-^+0(s) 
d\  oy 


Dropping  the  primes  to  simplify  the  notation: 


^Ix  +<t>Ixt  +s(^^2x  +<l>2xt  +<t>lx<t>lxx  j  =0(8^  )  (3-la) 

hly  +£(^2y  +^2yt  +<t>lx<l>lxy  =0(S^)  (3-lb) 


nu  +h(|)„,  +8(ti2,  +h<t)2,, +h({»,yy +(ri,(t)„X  +<t>ixhx)=  O(s')  (3-lc) 

where  subscripts  t,  x,  y,  x  indicate  derivatives. 

The  lowest  order  wave  field  (q,,<t),)  is  assumed  to  be  a  linear  superposition  of 
nearly  plane,  shoreward  propagating  waves  with  frequency  Qp  and  alongshore 
wavenumber 
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(3-2) 


00  ao  I 

S  7ap,,(x)exp[i(v|/p^(x)+-f^y-(Dpt)] 

p  =  -aD  q  =  -oo  ^ 

00  00  J 

=  Z  Z  ^ap,^(x)exp[i(v(/p^(x)+^<,y-capt)] 

p=-oo  q=-oo  p 

where  Op  =  pA©  and  =  qA^,  with  A©  and  A^  the  separation  of  adjacent  bands  in  the 
Fourier  representation,  and  the  spatial  variation  in  phase  is  given  by 


d\ 


ylW) 


+  sTp,qW 


(3-3) 


The  amplitude  (Up^ )  and  phase  (Vp,,)  functions  include  a  slow  variation  with  x,  owing 

to  shoaling,  refraction  and  non-linear  effects.  This  same  slow  space  variable  x  was  used 
previously  to  scale  the  beach  slope  in  Equations  (2-12)  and  (2-13).  In  order  for  r|,  and 
(j)i  to  be  real  requires: 


^p.q  ^-p.-q  ^  p.q  ^-p.-q  ^p.q  "^-p.-q 

Linear  terms  involving  q,  are  eliminated  from  the  governing  Equations  (3-la,b,c) 

by  cross  differentiating  (3-la)-—  (3-lc)and  —  (3-lb)  -  —  (3-lc)): 

ot  ox  dt  dy 

‘t>lxn  -(Hlxx)x  +s(^?2xt.  -(H2xx)x  +(‘|)lx<t>lxx),  "  ^  (h<|)  -('Hl't'lxX,  "  (hjCj),  J,)  =  0(8^  ) 

(3-4) 

<}>,>„  -(h(|),„),  -(h$2xx)y  +(<t>lx<l>,xy),  -^h^.xxyt.  -(hl<l>.xX,  -(hx‘t),x)y)  =0(8^) 

(3-5) 
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(3-7) 


where  F(t)  is  an  arbitrary  function  of  time  arising  from  the  integration  which  we  set  equal 
to  zero.  Substitution  of  the  lowest  order  wave  field  (Equation  (3-2)  and  (3-3))  into 
Equation  (3-7)  yields,  after  some  algebraic  manipulations: 


<l>2u-h«t'2xx  =  Z  S  4ap,,^+ap,,5Vh+i^a,,,VhTp,,--ap,,hco/  + 


'T^P.q - 

2  co„ 


1 1  ap-™,™-„a„.,„®pexpji(M/p_^,_„+M^„.„-v|;p,)]Uxpji(v);p^(x)+f^y-c0pt)]+O(E) 


The  left-hand  side  of  Equation  (3-8)  is  the  wave  equation  for  (t)2 .  To  prevent 

resonant  growth  of  <|)2  (which  would  upset  the  perturbation  expansion,  Equation  (2-7)) 
requires  that  the  forcing  terms  on  the  right-hand  side  of  Equation  (3-8)  do  not  contain  any 
free  oscillations.  Since  all  these  terms  obey  the  lowest-order  dispersion  relation 
(Equation  (3-3)),  it  follows  that  the  right-hand  side  of  Equation  (3-8)  must  vanish. 
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Collecting  like  frequencies  and  alongshore  wavenumbers  and  solving  for  the  real  and 
imaginary  parts  yields  a  coupled  set  of  equations  for  the  amplitude  and  phase  evolution  of 
the  lowest  order  waves: 


a  ~  =  — a 

p.qx  4 


3  -- 
-h 

8  ^ 


CO  CO 

^p-m.q-n^m.n  m.n  ~  H^p.q  ) 


(3-9a) 


=-Vh«; 


I  Z 


^p-m.q-n^m.n  p-m,q-n  m.n  p 


(3-9b) 


The  first  term  on  the  right  hand  side  of  Equation  (3-9a)  is  the  amplitude  growth  of 
the  mode  due  to  shoaling.  The  first  and  second  terms  on  the  right  hand  side  of  Equation 
(3 -9b)  are  dispersion  and  refraction  corrections  to  the  wave  phase.  The  double 
summation  terms  give  the  amplitude  growth  and  phase  changes  of  a  mode  due  to  resonant 
non-linear  interactions  of  all  possible  triads  in  which  the  mode  participates.  For  uni¬ 
directional,  normally  incident  waves,  Equations  (3-9a,b)  reduce  to  Freilich  and  Guza’s 
(1984)  “consistent  shoaling  model”  (Equations  (16a,b)  in  their  paper). 

Finally,  we  re-dimensionalize  Equations  (3-9a,b),  and  re-express  them  in  unsealed 
(physical)  coordinates.  Dropping  the  used  previously  to  denote  dimensional 
coordinates,  and  using  Equation  (3-3),  we  obtain 
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‘ ''"8»,., "  ’  7J.?.  2  *'^"--Vm) 


(3-lOb) 
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IV.  COMPARISON  TO  FINITE  DEPTH  THEORY 

In  this  Chapter  the  amplitude  and  phase  evolution  given  by  Equation  (3-10a,b)  are 
compared  to  predictions  of  second-order  nonlinear  finite  depth  theory  (Hasselmarm, 
1962).  The  Boussinesq  equations  were  integrated  over  a  plane  beach  for  a  single  resonant 
triad  consisting  of  two  primary  waves  incident  from  deep  water  that  drive  a  secondary 
wave  component.  Both  sum-  and  difference-interactions  are  compared  for  a  range  of 
commonly  observed  beach  slopes,  incident  wave  amplitudes,  frequencies,  and 
propagation  directions.  The  objective  of  these  comparisons  is  to  examine  the  similarities 
and  differences  of  Boussinesq  and  finite  depth  theory  predictions  in  the  transition  region 
from  deep  to  shallow  water,  and  in  particular  to  investigate  the  dependence  of  the 
secondary  wave  growth  on  the  directional  spreading  angle  of  the  primary  waves.  While 
this  approach  is  useful  to  describe  the  initial  growth  of  a  secondary  wave  component,  in 
very  shallow  water  the  secondary  wave  amplitude  grows  to  become  comparable  to  the 
primary  wave  amplitudes  and  thus  fiilly  coupled  spectral  calculations  are  needed  to 
predict  the  nonlinear  shoaling  evolution. 

For  a  single  triad,  we  can  simplify  the  notation  used  in  Equation  (3-10a,b)  by 
introducing  a  single  index  i=l,2,3,  where  mode  i  has  a  frequency  co^'^ 
(0  <0)*’’  alongshore  wavemunber  ,  amplitude  a^'^ ,  and  phase  v|/^'^ ,  and 

the  triad  satisfies  the  resonance  rule: 
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The  amplitude  and  phase  evolution  equations  for  this  triad  can  be  expressed  (neglecting 
energy  exchanges  with  any  other  components)  compactly  as: 

a<‘>  =S^’’+N“*sinB“^  (4-la) 


(shoaling  term) 


(shallow  water  wave  number) 


(dispersion  term) 


(refraction  term) 


(nonlinear  interaction  term) 


(biphase) 


A  first-order  finite  difference  scheme  using  an  Arakawa  A  grid  with  a  step  size  of 
0.1  m  was  used  to  integrate  Equation  (4-1).  Test  runs  with  grid  spacings  of  0.01  m  and 
1.0  m  yielded  identical  results,  confirming  the  accuracy  of  the  numerical  scheme. 
Additional  tests  of  the  numerical  scheme  include  a  verification  of  energy  conservation. 
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In  second-order  finite  depth  theory,  two  primary  waves  with  frequencies  and 
(0^^^  (co^'^  >  (0^^^  >  0)  and  alongshore  wavenumbers  t''’  and  excite  secondary  wave 
components  with  the  sum-frequency  and  wavenumber  (co^'^  +  and  the 

difference-frequency  and  wave  number  ).  In  finite  depth  theory,  the 

free  surface  elevation  is  given  by 

q(x,y,t)=a^'^cos(v|/^'’(x)+^‘'V-®^'’t)+a‘^^cos(v|/^^^(x)+^‘^V-®^^^t) 
+C^a‘'*(x)a‘^*(x)cos^^'*(x)+n/^^^(x)+|j?‘'^+^‘^’ Jy-[co*'^-i-co^^’ 
+C.a^'’(x)a‘^’(x)cos^^''(x)-v|/‘^’(x)+[^^‘'’-^‘^’ Jy-[co‘'*-co*^* 

+ other  second -order  terms 

where  the  primary  wave  amplitudes  a^'^,  a^^^  and  phases  are  given  by  the  linear 

shoaling  and  refraction  relations: 


a'‘’(x)  = 


Cj'j^cosG^ 


(i) 


Cg’(x)  cos(0‘‘'(x)^ 


.(i) 


M/<‘’(x)  =  v|;<‘>  + Jk<'>(Qd;  ;  i  =  l,2  (4-3) 


with  subscripts  o  indicating  deep  water  values  at  x  =  0.  Cg  denotes  the  group  velocity  and 
0^''  is  given  by  Snell’s  law.  The  coupling  coefficients  C+  and  C.  for  the  sum-  and 
difference-interactions  are  given  by  (Hasselmann,  1962): 
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gk*'^  tanh^k*'^  ±k‘'’[h^^(o*'' 


gk<').k(2) 

g 

k<'> 

|2  1  ,-,x|2  A) 

k® 

■¥  ' 

g 

2(cd<'>±co'^') 

^^co^'^cosh^ 

(jj(l)^  ±  (lj(l)gj(2)  ^  ^(2)^  gl^(l)  _  1^(2) 

2g  2ol)''’co‘^* 

with  k^'^  the  vector  wavenumber  of  component  i  (k*‘\x),  i 

Finite  depth  theory  predictions  of  primary  and  secondary  wave  amplitudes  in  10 
m  depth  were  used  to  initialize  the  Boussinesq  model  (Equation  (4-1)).  The  biphase  B  in 
Equation  (4-1)  was  matched  to  the  finite  depth  value  (0  or  tt,  depending  on  the  sign  of  C+, 
Equations  (4-2),  (4-4)). 

Boussinesq  and  finite  depth  theory  predictions  of  sum-  and  difference-frequency 
waves  forced  by  a  pair  of  swell  components  with  frequencies  0.08  and  0.12  Hz  on  a  slope 
0.01  beach  are  compared  in  Figures  2  and  3.  The  deep  water  incident  wave  amplitudes  in 
these  simulations  are  0.2  m.  In  Figure  2  results  are  shown  for  two  normally  incident 
swell  components.  In  Figure  3  the  0.08  Hz  component  is  normally  incident  while  the 
0. 12  Hz  component  has  a  deep  water  incidence  angle  of  60°.  In  all  cases  the  Boussinesq 
and  finite  depth  theory  predictions  show  very  similar  secondary  wave  growth  between  10 
and  4  m  depth.  Overall,  the  amplitudes  predicted  by  finite  depth  theory  are  slightly 
higher  than  those  predicted  by  Boussinesq  theory,  suggesting  that  finite  depth  theory 
overpredicts  the  secondary  wave  growth  when  the  triad  is  close  to  resonance.  Boussinesq 
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theory  predictions  of  the  biphase  are  generally  close  to  the  finite  depth  theory  values  of  0 
(sum-interactions)  and  n  (difference-interactions).  The  secondary  wave  amplitude 
growth  in  the  Boussinesq  predictions  is  controlled  by  the  biphase  of  the  triad  (Equation 
(4-12)).  The  growth  vanishes  when  B  is  equal  to  the  finite  depth  values  of  0  or  n.  Hence, 
the  Boussinesq  and  finite  depth  theories  do  not  smoothly  match  in  10  m  depth  and  the 
Boussinesq  biphase  appears  to  oscillate  around,  and  converge  to  the  small  value  which 
causes  the  secondary  wave  to  grow  at  the  same  rate  as  predicted  by  finite  depth  theory. 

The  amplitude  and  biphase  discrepancies  between  Boussinesq  and  finite  depth 
theory  predictions  are  generally  larger  for  difference-interactions  than  for  sum- 
interactions  (Figures  2,  3  and  other  examples  discussed  below).  Errors  are  expected  in 
both  finite  depth  and  Boussinesq  theory  predictions  for  triads  involving  a  low 
(infragravity)  frequency  secondary  wave  component  because  the  infragravity  wavelengths 
are  much  longer  than  those  of  the  incident  swell.  For  example,  in  the  model  run  shown  in 
Figure  (2a,c)  the  wavelength  of  the  0.04  Hz  difference-frequency  component  varies 
between  250  and  150  m  and  thus  the  depth  variations  over  a  wavelength  (1.5  to  2.5  m) 
are  appreciable. 

While  the  predicted  sum-frequency  secondary  wave  amplitudes  are  only  slightly 
smaller  for  obliquely  propagating  primary  waves  than  for  normally  incident  primary 
waves  (cf  Figures  2b  and  3b),  the  60°  angular  separation  of  the  primary  waves  causes 
about  an  order  of  magnitude  reduction  in  the  difference-frequency  secondary  wave 
amplitude  (cf  Figures  2a  and  3a).  The  dependence  of  the  sum-  and  difference-frequency 
secondary  wave  amplitudes  on  primary  wave  incidence  angles  is  illustrated  in  Figures  4- 
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6.  Secondary  wave  amplitudes  in  4  m  depth,  normalized  by  the  product  of  the  primary 
wave  amplitudes  (in  the  same  depth),  are  shown  versus  the  directional  spreading  angle  of 
the  primary  waves  in  deep  water  for  a  range  of  beach  slopes  and  primary  wave 
frequencies  and  amplitudes.  In  finite  depth  theory  the  normalized  amplitude  is  equal  to 
the  coupling  coefficient  C  (Equation  (4-2))  and  thus  can  be  interpreted  as  the  secondary 
wave  response. 

In  Figure  4  results  are  shown  for  primary  wave  pairs  with  frequencies  (0.09  Hz, 
0.1 1  Hz)  and  (0.07  Hz,  0.13  Hz),  and  amplitudes  0.2  m  on  a  slope  0.01  beach.  In  Figure 
5  the  response  of  a  (0.08  Hz,  0.12  Hz)  primary  wave  pair  with  amplitudes  0.1  m  is 
compared  to  the  response  of  the  same  frequency  pair,  but  with  larger  amplitudes  0.4  m, 
on  a  slope  0.01  beach.  Results  for  0.2  m  amplitude  primary  waves  with  frequencies  (0.08 
Hz,  0.12  Hz)  on  two  different  beach  slopes  0.003  and  0.03  are  shown  in  Figure  6.  The 
Boussinesq  and  finite  depth  theory  predictions  of  the  dependence  of  the  secondary  wave 
response  on  the  directional  spreading  angle  of  the  two  primary  waves  are  generally  in 
good  agreement.  In  all  cases  sum-interactions  are  insensitive  to  the  spreading  angle  with 
approximately  a  10-20%  decrease  in  the  response  when  the  spreading  angle  increases 
from  0  to  60°.  The  magnitude  of  the  predicted  sum-frequency  response  is  insensitive  to 
the  primary  wave  frequencies  and  the  beach  slope  with  Boussinesq  predictions 
consistently  about  25%  lower  than  finite  depth  predictions  (cf.  Figures  4b,d;  6b,d). 
However,  while  Boussinesq  predictions  of  the  sum-frequency  response  for  0.1  m 
amplitude  primary  waves  (Figure  5b),  are  in  very  close  agreement  (about  10-20%  lower) 
with  the  finite  depth  values,  the  Boussinesq  predictions  for  0.4  m  amplitude  waves  are 
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about  30-50%  lower  than  the  finite  depth  values  (Figure  5d).  In  this  case  finite  depth 
theory  predicts  secondary  wave  amplitudes  that  are  comparable  to  primary  wave 
amplitudes  and  thus  is  not  valid.  Boussinesq  theory  predicts  a  reduced  response  for  these 
large  amplitude  primary  waves  that  is  virtually  independent  of  the  spreading  angle. 

The  difference-frequency  secondary  wave  response  is  much  more  sensitive  to  the 
directional  spreading  angle  than  the  sum-frequency  response  (cf  Figures  4a, c;  5a, c;  6a,c 
to  Figures  4b, d;  5b,d;  6b,d).  An  increase  in  deep  water  spreading  angle  from  0  to  60° 
causes  a  reduction  in  the  difference-frequency  response  by  a  factor  3  to  10.  The 
Boussinesq  predictions  are  generally  in  good  agreement  with  finite  depth  theory 
predictions,  showing  somewhat  lower  response  levels  but  approximately  the  same 
variation  with  directional  spreading  angle.  However,  the  agreement  is  much  better  on  a 
gentle  (0.003)  slope  (Figure  6a)  than  on  a  steep  (0.03)  slope  (Figure  6c).  In  the  latter  case 
the  depth  changes  from  10  to  4  m  over  a  distance  comparable  to  the  secondary 
wavelength  and  the  small  beach  slope  assumption  used  in  both  Boussinesq  and  finite 
depth  theory  is  grossly  violated.  Similar,  but  less  pronounced  discrepancies,  between  the 
Boussinesq  and  finite  depth  theory  predictions  for  0.02  Hz  secondary  waves  (Figure  4a, 
with  wavelengths  0(500  m))  suggests  a  break  down  of  the  small  beach  slope 
approximation. 

Overall,  the  comparisons  of  Boussinesq  and  finite  depth  theory  predictions  of 
secondary  wave  growth  in  the  transition  region  from  deep  to  shallow  water  show 
surprisingly  good  agreement.  The  Boussinesq  predictions  of  both  sum-  and  difference- 
frequency  secondary  wave  amplitudes  are  slightly  lower  than  finite  depth  predictions  but 
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the  dependence  on  directional  spreading  angle  predicted  by  both  theories  are  in  excellent 
agreement. 
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V.  CONCLUSIONS 


Although  wave  transformation  on  beaches  is  generally  well  described  by  one¬ 
dimensional  models,  the  directionality  of  waves  is  of  crucial  importance  to  a  variety  of 
nearshore  processes  including  inffagravity  waves,  longshore  currents  and  sediment 
transport.  In  this  paper  Freilich  and  Guza’s  (1984)  one-dimensional  spectral  model  for 
the  shoaling  of  imi-directional  waves  on  a  beach,  based  on  the  Boussinesq  equations  for 
weakly  nonlinear,  weakly  dispersive  waves  in  slowly  varying  depth  (Peregrine,  1 967)  is 
extended  to  two  dimensions.  A  beach  with  straight  and  parallel  depth  contours  is 
assumed  on  which  wave  incidence  angles  are  reduced  by  refraction.  In  this 
approximation,  dispersion  and  directional  spreading  effects  are  shown  to  be  of  the  same 
order,  and  near-resonant  nonlinear  interactions  occur  between  any  pair  of  incident  wave 
components.  A  two-dimensional  Fourier  representation  in  the  frequency  and  alongshore 
wavenumber  domain  is  used  to  describe  the  lowest  order  wave  field,  yielding  a  coupled 
set  of  equations  for  the  amplitude  and  phase  evolution  of  each  Fourier  mode.  For  uni¬ 
directional  waves  these  equations  reduce  to  Freilich  and  Guza’s  model. 

Model  results  for  the  simple  case  of  a  single  triad  of  waves  are  compared  to 
second-order  nonlinear,  fully  dispersive  finite  depth  theory  (Hasselmann,  1962)  for  a 
typical  range  of  beach  slopes  and  swell  amplitudes,  frequencies  and  propagation 
directions.  In  these  model  runs  two  primary  swell  components  force  a  third  (secondary) 
wave  component  with  the  sum-  or  difference-frequency  and  alongshore  wavenumber. 
The  Boussinesq  model  (not  valid  in  deep  water)  was  initialized  using  the  finite  depth 
theory  predictions  of  the  amplitudes  and  biphase  in  10  m  depth.  Boussinesq  and  finite 
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depth  theory  predictions  of  both  sum-  and  difference-interations  are  compared  in  the 
transition  region  from  deep  to  shallow  water  (10-4  m  depth)  where  both  theories  are 
expected  to  be  reasonably  accurate. 

In  general,  both  theories  predict  very  similar  secondary  wave  amplitude  growth, 
with  Boussinesq  predictions  typically  slightly  lower  than  finite  depth  predictions.  The 
dependence  of  the  secondary  wave  response  on  the  directional  spreading  angle  of  the 
incident  swell  components  in  deep  water  predicted  by  both  theories  are  in  excellent 
agreement.  Whereas  the  sum-frequency  response  is  only  slightly  reduced  for  large 
spreading  angles,  difference-interactions  are  sensitive  to  the  primary  wave  directions  with 
typical  reductions  in  secondary  wave  amplitudes  of  a  factor  of  3  to  10  when  the  primary 
wave  spreading  angle  increases  from  0°  to  60°. 

Not  surprising,  discrepancies  between  Boussinesq  and  finite  depth  theory  are  most 
pronounced  for  low  (infragravity)  frequency  secondary  waves  on  a  relatively  steep  beach. 
In  this  case  depth  variations  are  appreciable  over  distances  comparable  to  the  secondary 
wavelength  and  the  small  beach  slope  assumption  used  in  both  theories  is  violated. 
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APPENDIX 


Figure  1.  Definitional  sketch  of  variables  and  coordinate  frame.  Obliquely 
propagating  waves  shoal  upon  a  beach  having  straight  and  parallel  depth  contours,  ho  is  a 
representative  water  depth  in  the  region,  0  is  the  wave  propagation  angle,  and  k  is  the 
wavenumber  vector. 

Figure  2.  Comparisons  of  Boussinesq  and  finite  depth  theory  predictions  of  the 
shoaling  evolution  of  a  single  wave  triad.  A  pair  of  normally  incident  primary  swell 
components  with  frequencies  0.08  Hz  and  0.12  Hz  drives  a  difference-frequency  (0.04 
Hz)  secondary  wave  (left  panels)  and  a  sum-frequency  (0.20  Hz)  secondary  wave  (right 
panels).  The  upper  and  lower  panels  show  the  evolution  of  the  amplitudes  and  biphase. 
Deep  water  primary  wave  amplitudes  are  0.2  m. 

Figure  3.  Same  model  comparisons  as  shown  in  Figure  2  but  with  an  obliquely 
propagating  0.12  Hz  primary  wave  component  (deep  water  incidence  angle  of  60°). 

Figure  4.  Secondary  wave  response  (defined  as  the  secondary  wave  amplitude 
normalized  by  the  product  of  the  primary  wave  amplitudes)  in  4  m  depth  versus  the  deep 
water  directional  spreading  angle  of  the  primary  wave  components.  The  primary  wave 
frequencies  are  0.09  Hz,  0.11  Hz  (upper  panels)  and  0.07  Hz,  0.13  Hz  (lower  panels). 
The  lower  frequency  primary  wave  component  is  normally  incident  while  the  deep  water 
incidence  angle  of  the  higher  frequency  primary  wave  varies  between  0  and  60°. 
Difference  and  sum  interactions  are  shown  in  the  left  and  right  panels,  respectively.  The 
beach  slope  is  0.01  and  primary  wave  amplitudes  in  deep  water  are  0.2  m. 
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Figure  5.  Secondary  wave  response  in  4  m  depth  for  primary  waves  with 
frequencies  0.08  and  0.12  Hz  on  a  beach  with  slope  0.01.  Upper  and  lower  panels  show 
results  for  deep  water  primary  wave  amplitudes  of  0.1  and  0.4  m  respectively  (same 
format  as  Figure  4). 

Figure  6.  Secondary  wave  response  in  4  m  depth  for  0.2  m  amplitude  primary 
waves  with  frequencies  0.08  and  0.12  Hz  on  beaches  with  slope  0.003  (upper  panels)  and 
0.03  (lower  panels)  (same  format  as  Figure  4). 


32 


BIPHASE  (deg)  AMPLITUDE  (m) 


BIPHASE  (deg)  AMPLITUDE  (m) 


10  8  6  4 


DEPTH  (m) 

Igure  3 


35 


DIFFERENCE  INTERACTIONS 


SUM  INTERACTIONS 


DIFFERENCE  INTERACTIONS 


SUM  INTERACTIONS 


LIST  OF  REFERENCES 


Beach,  R.,  and  Sternberg,  R.,  1988,  “Suspended  Sediment  in  the  Surf  Zone:  Response  to 
Cross-shore  Infragravity  Motion”,  Marine  Geology  80,  61-79. 

Boussinesq,  J.,  1871,  Theorie  Generate  des  Mouvements  qui  Sont  Propages  dans  un 
Canal  Rectangulaire  Horizontal,  C.R.  Hebd.  Seane,  Acad.  Sci.,  Paris  73, 256- 
260. 

Elgar,  S.,  Herbers,  T.H.C.,  and  Guza,  R.T.,  1994  “Reflection  of  Ocean  Surface  Gravity 
Waves  from  a  Natural  Beach”,  J.  Phys.  Oceanogr.,  24(7),  1503-1511. 

Elgar,  S.  and  Guza,  R.T.,  1985,  “Observations  of  Bispectra  of  Shoaling  Surface  Gravity 
Waves,”  J.  Fluid  Meek,  161,425-448. 

Elgar,  S.  and  Guza,  R.T.,  1986,  “Nonlinear  Model  Predictions  of  Bispectra  of  Shoaling 
Surface  Gravity  Waves”,  J.  Fluid  Meek,  167,  1-18. 

Freilich,  M.H.,  and  Guza,  R.T.,  1984,  “Nonlinear  Effects  on  Shoaling  Surface  Gravity 
Waves,”  Phil.  Trans.  R.  Soc.  Land,  A31 1,  1-42. 

Hasselmann,  K.,  1 962,  “On  the  Nonlinear  Energy  Transfer  in  a  Gravity-Wave  Spectrum, 
1,  General  Theory,”  J.  Fluid  Meek,  12, 481-500. 

Herbers,  T.H.C.,  and  Guza,  R.T.,  1991,  “Wind-Wave  Nonlinearity  Observed  at  the  Sea 
Floor,  Part  1,  Forced  Wave  Energy,”  J.  Phys.  Oeeanogr.,  21(12),  1740-1761. 

_ and _ ,  1 992,  “Wind-Wave  Nonlinearity  Observed  at  the  Sea 

Floor,  Part  2,  Wavenumbers  and  Third-Order  Statistics,”  J.  Phys.  Oeeanogr., 
22(5),  489-504. 

_ and _ ,  1994,  “Nonlinear  Wave  Interactions  and  High-Frequency 

Seafloor  Pressure”,  J.  Geophys.  Res.,  99(c5),  10035-10048. 

Herbers,  T.H.C.,  Lowe,  R.L.,  and  Guza,  R.T.,  1992,  “Field  Observations  of  Orbital 

Velocities  and  Pressure  in  Weakly  Nonlinear  Surface  Gravity  Waves,”  J.  Fluid 
Meek,  245,413-425. 

Herbers,  T.H.C.,  Elgar,  S.,  and  Guza,  R.T.,  1994,  “Infragravity-Frequency  (0.005-0.05 

Hz)  Motions  on  the  Shelf,  Part  I,  Forced  Waves,”  J.  Phys.  Oeeanogr.,  24(5),  917- 
927. 


39 


Berbers,  T.H.C.,  Elgar,  S.,  Guza,  R.T.,  and  O’Reilly,  W.C.,  1995a,  “Infragravity- 

Frequency  (0.005-0.05  Hz)  Motions  on  the  Shelf,  Part  II;  Free  Waves”,  J.  Phys. 
Oceanogr.,  in  press. 


Herbers,  T.H.C.,  Elgar,  S.,  and  Guza,  R.T.,  1995b,  “Generation  and  Propagation  of  Surf 
Beat”,  J.  Geophys.  Res.,  submitted. 

Holman,  R.A.,  Huntley,  D.A.,  and  Bowen,  A.J.,  1978,  “Infragravity  Waves  in  Storm 

Conditions,”  Proc.  16th  Conf.  on  Coastal  Engineering,  American  Society  of  Civil 
Engineers,  268-284. 


Liu,  P.L.F.,  Yoon,  S.B.,  and  Kirby,  J.T.,  1985,  “Nonlinear  Refraction-Diffraction  of 
Waves  in  Shallow  Water,”  j;  Fluid  Meek,  153, 185-201. 

Longuet-Higgins,  M.S.  and  Stewart,  R.W.,  1962,  “Radiation  Stress  and  Mass  Transport 
in  Gravity  Waves,  with  Application  to  ‘Surf  Beats’,”  7!  Fluid  Meek,  13, 481-504. 

Mei,  C.C.  and  LeMehaute,  B.,  1966,  “Note  on  the  Equations  of  Long  Waves  over  an 
Uneven  Bottom,”  J.  Geophys.  Res.,  91(B7),  7343-7358. 

Munk,  W.H.,  1949,  “Surf  Beats,”  Trans.  Amer.  Geophys.  Un.,  30,  849-854. 

Okihiro,  M.,  Guza,  R.T.,  and  Seymour,  R.J.,  “Bound  Infragravity  Waves”,  J.  Geophys. 
Res.,  91,  11,453-11,469,  1992. 

Peregrine,  D.H.,  1967,  “Long  Waves  on  a  Beach,”  J.  Fluid  Meek,  27,  Part  4,  815-827. 

Phillips,  O.M.,  1960,  “On  the  Dynamics  of  Unsteady  Gravity  Waves  of  finite  Amplitude, 
Part  1,  The  Elementary  Interactions,”  J.  Fluid  Meek,  9,  193-217. 

Stokes,  G.G.,  1 847,  “On  the  Theory  of  Oscillatory  Waves.”  Trans.  Camb.  Phil.  Soe  ,  8, 
441-455. 


Tucker,  M.J.,  1950,  “Surf  Beats:  Sea  Waves  of  1  to  5  Minutes  Period,”  Proe.  Roy  Soe 
A202,  565-573. 

Ursell,  F.,  1953,  “The  Long-Wave  Paradox  in  the  Theory  of  Gravity  Waves,”  Proe. 
Camb.  Phil.  Soe.,  49,  685-694. 

Wright,  L.D.,  Guza,  R.T.,  and  Short,  A.D.,  1982,  “Dynamics  of  a  High-Energy 
Dissipative  Surf  Zone,”  Mar.  GeoL,  45,  41-62. 


40 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center . 2 

Cameron  Station 

Alexandria,  VA  22304  -  6145 

2.  Library,  Code  52 . 2 

Naval  Postgraduate  School 

Monterey  CA  93943  -  5101 

3.  Professor  T.H.C.  Berbers,  Code  OC/He . 8 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey  CA  93943  -5121 

4.  Professor  Frenzen,  Code  MA/Fr . 1 

Department  of  Mathematics 

Naval  Postgraduate  School 
Monterey  CA  93943  -  5121 

5.  Professor  E.B.  Thornton,  Code  OC/Tm . 1 

Naval  Postgraduate  School 

Monterey  CA  93943  -  5121 

6.  Professor  T.C.  Lippmann,  Code  OC/Li . 1 

Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey  CA  93943  -  5121 

7.  Mark  Burton . 2 

1590  Harbor  Blvd. 

Belmont,  CA  94002 


41 


